Depairing currents in superconducting films of Nb and amorphous MoGe. 
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We report on measuring the depairing current Jdp in thin superconducting films as a function of 
temperature. The main difficulties in such measurements are that heating has to be avoided, either 
due to contacts, or to vortex flow. The latter is almost unavoidable since the sample cross-section is 
usually larger than the superconducting coherence length and the magnetic field penetration depth 
As. On the other hand, vortex flow is helpful since it homogenizes the distribution of the current 
across the sample. We used a pulsed current method, which allows to overcome the difficulties 
caused by dissipation and measured the depairing current in films of thin poly crystalline Nb (low 
As, low specific resistance p) and amorphous Moo.tGco.s (high As, high p), structured in the shape 
of bridges of various width. The experimental values of Jdp for different bridge dimensions are 
compared with theoretical predictions by Kupriyanov and Lukichev for dirty limit superconductors. 
For the smallest samples we find a very good agreement with theory, over essentially the whole 
temperature interval below the superconducting critical temperature. 
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I. INTRODUCTION 

The superconducting current density Jg is a unique 
feature of a superconducting material. It can be ex- 
pressed as Js=engVg, where rig and Vg are the density 
and velocity of the superconducting electrons respec- 
tively, and e is the electron charge. Increasing 3g leads 
to increase of Vg but also to a reduction of the number 
of Cooper pairs. Finally, when Jg reaches the depair- 
ing current Jdp, the amount of carriers is not enough to 
support the supercurrent and the superconducting state 
collapses. For conventional superconductors the temper- 
ature dependence of Jdp near the critical temperature 
Tc is given by the classical Ginzburg-Landau (GL) ex- 
pression Jg^{t) = J|,^(0)(1 - <)3/2, where t = T/T^, 
and Jdj^{0) is the depairing current at zero temperature. 
Early work on determining Jdp in Sn microbridges can 
be found in P, Q . The GL approach becomes invalid at 
lower temperatures, since the conditions 1 — T/Tc 

for clean limit superconductors {k is Ginzburg-Landau 
parameter), or (Tc — T)<^ Tc for dirty limit supercon- 
ductors, are no longer fulfilled. A more complete and 
quantitative theory, valid for all temperatures and arbi- 
trary mean free path, was developed by Kupriyanov and 
Lukichev (KL), who obtained the numerical solution of 
the Eilenberger equations for a superconductor carrying 
a current, with the velocity of the Cooper-pairs propor- 
tional to a phase gradient of the superconducting order 
parameter A Notably, their theory gives the same 
expression for Jdp(t) as GL theory for the temperature 
region close to Tc and also yields the correct expressions 
for Jdp(O) in terms of the materials constants. 
The amount of theoretical work done on depairing cur- 
rents in conventional superconductors contrasts sharply 
with a lack of experimental observations. A major is- 
sue here is the requirement with respect to sample di- 
mensions. In principle, the sample width should not be 
larger than both the penetration depth As, and the co- 
herence length ^s- The first condition is needed to avoid 
current piling up at the edges, because of the Meiss- 



ner effect 0. For a superconducting film A^ is given by 
Xl/dg, (dg <C Afc) where A^ is the bulk London penetra- 
tion depth, dg is film thickness, and the magnetic field 
is taken perpendicular to the film plane. At low tem- 
peratures in case of dirty superconductors it becomes 
Xl{£,o/£ds), where is the BCS coherence length, and 
£ is the elastic mean free path. A typical value of Xb 
for instance for polycrystalline Nb is 50 nm; for amor- 
phous materials such as a-Moo.rCeo.s, which will also 
be discussed below, Af, is much larger, of the order of 
0.5 /im. The condition on ^g must be fulfilled when vor- 
tex nucleation and flow is to be prevented, which cause 
dissipation in sample before the Jdp is reached. Exact 
calculations made by Likharev 5] show that the smallest 
sample width below which no vortex can appear equals 
4.4^3 (T), where ^s{T) is the Ginzburg-Landau coherence 
length given by Cs{T) = 0.85 Cs(0)/VT^, with ^,(0) = 
-v/^qI. Typical values of ^s(O) for our Nb and Moo.yGeo.a 
are 12 nm (because of the small mean free path) and 
5 nm, respectively. The only case where both conditions 
can be implemented is a thin aluminum film shaped in 
a form of a narrow (about l/im) bridge. The BCS co- 
herence length for Al is of the order of 1.5 /im, while the 
penetration depth can be increased to a similar value by 
decreasing the film thickness. Romijn et al. @ showed 
that for such system the experimental values of the de- 
pairing current density were in excellent agreement with 
KL theory for temperatures down to 0.2t. In case of Nb 
and Moo.rGeo.s films one would have to go to a bridge 
width not larger than 30 nm in order to prevent vortex 
appearance. 

However, vortex motion also has an advantage, since it 
will homogenize the current distribution 7]. The main 
problem then in determining Jdp is to avoid sample heat- 
ing, cither by dissipation due to vortex motion, or e.g. 
to heating in the contacts due to the relatively large cur- 
rents. In this paper we demonstrate that the undesired 
sample heating can be avoided by using a pulsed current 
method. We use different superconductors, with widely 
different values of Jdp- Specifically, we use Nb with low 
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Xb and also relatively low specific resistance p (around 
7 /iOcm) and amorphous (a-)Moo.7Geo 3 with large 
and a large p « 160 fiQcm. Especially the large p easily 
leads to dissipation in the neighborhood of the transi- 
tion to the normal metal state. Films of different thick- 
nesses were patterned into bridges of different width w^. 
The experimental values we obtain for the depairing cur- 
rent density Jdp{t) are in very good agreement with the 
KL calculations, assuming that the current distribution 
across the samples is perfectly homogeneous. 



II. EXPERIMENT 

Nb single layer films were grown by dc magnetron 
sputtering in an ultra high vacuum system with a back- 
ground pressure of about 10~^° mbar and an Ar sput- 
tering pressure of 6x10"'^ mbar. Films of a-Moo.yGeo.a 
were deposited in a RF-diode sputtering system with a 
background pressure of 10^^ mbar in an Ar pressure of 
8x10""^ mbar. Sputtering rates for Nb and a-Moo.rGeo.a 
were 0.8 A/s and 1.2 A/s respectively. Both materi- 
als were grown on Si (100) substrates. The thickness 
of the films was determined during the deposition by a 
crystal thickness monitor, which was calibrated by low 
angle X-ray diffraction measurements and Rutherford 
Backscattering. For the depairing current experiments, 
samples were structured in the shape of strips of differ- 
ent cross-section by e-beam lithography and Ar-ion etch- 
ing. The structure included the contacts. In the case of 
a-Moo.yGeo.s, samples were water-cooled during deposi- 
tion and liquid nitrogen-cooled during etching, in order 
to prevent undesirable film crystallization. The typical 
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FIG. 1: Sample layout. The measurement procedure was per- 
formed with a classical 4-point scheme. The massive current 
leads provide a good heat sink. 

geometry of the samples is shown in Fig. 1. In all cases 



the distance between voltage leads was 100±1 pm. The 
width of resistive transition from the normal into the su- 
perconducting state was about 30 mK for all samples. 
An example for both materials is given in Fig. 2. Trans- 
port measurements in the normal state yielded an aver- 
age value of specific resistance p of about 160 pflcm for 
MoQ.rGeo.s and 7.2 pflcm for Nb samples respectively. 
For a-Moo.7Geo.3 the elastic mean free path i is taken to 
be 0.4 nm |^, of the order of the interatomic distances 
and these samples are clearly in the dirty limit. For Nb, 
using the expressions of the free electron model with the 
product p£ = 3.75 xlO~^^ i7m^ and the Fermi velocity 
vp = 5.6 xlO^ m/s we find £ = 5.2 nm. Comparing 
this value to 1^0= 39 nm for Nb it is seen that the 
dirty limit condition ^ ^ ^0 is also satisfied. The de- 
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FIG. 2: Resistance normalized to its normal state value at 10 
K as a function of temperature for a Nb bridge (wa — 1 ^m, 
ds = 20 nm) and an a-Moo.rGeo.s bridge (ws = 2 pm, ds=64 
nm) 

pairing currents measurements were performed in a ^He 
cryostat shielded from external magnetic fields by a long 
permalloy (Nio.8Feo.2) screen annealed in hydrogen at- 
mosphere. Hall probe measurements showed a constant 
magnetic field background less than 10~^ T. The samples 
were mounted on a massive brass holder with a resistive 
heater. In order to reduce possible errors in the temper- 
ature determination because of the temperature gradient 
along the sample holder, all samples were placed in im- 
mediate proximity to the thermometer. The temperature 
stability during the experiment was about 1 mK . For 
determination of the critical current value Idp at differ- 
ent temperatures a pulsed current method was used, in 
which current pulses with a growing amplitude were sent 
through the sample. The average duration of a single 
pulse was about 3.00±0.05 ms. Each pulse was followed 
by a long pause of 7.0±0.1 s. The voltage response of the 
system was observed on an oscilloscope triggered for the 
time of a single pulse. To improve the signal resolution a 
differential amplifier was used, combined with low-noise 
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band filters. A typical current(/) - voltage(y) charac- 
teristic for a-Moo.rGeo.s at a reduced temperature oft — 
0.74 is shown in Fig. 3. One can see a clear jump from the 
superconducting to the normal state at Idp- For temper- 
atures close to Tc a small onset of voltage was observed 
in all samples, probably because of vortex motion. In 
order to make certain that this effect has no influence on 
the determination oi Idp, the temperature was monitored 
during every current pulse. Measurable differences were 
found very close to Idp, as shown in Fig. 3. We conclude 
that a short pulse in a combination with a long pause 
does not cause sample heating and keeps the system in 
temperature equilibrium until the dissipation related to 
the normal state occurs. 
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FIG. 3: Typical dependence of voltage V (open circles) and 
temperature T (open stars) on current /, measured on a 2 jj,m 
wide a-Moo.rGeo.s bridge. 
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FIG. 4: Experimental results for pair-braking current Jdp as 
function of reduced temperature for a Nb bridge (ds=20 nm, 
Ws = l/im) and an a-Moo.rGeo.s bridge (ds=64 nm, W3=2/im). 



Sample 


d. 
[nm] 




[K] 


P 

[^r2*cm] 


Jdp(O) 
10"[A/m2] 


jfp^(O) 
10"[A/m2] 


Nb 


20 


1.0 


8.3 


7.25 


17 


15 


Nb 


40 


2.0 


9.0 


7.24 


16 


17 


Nb 


53 


2.5 


9.0 


7.24 


19 


17 


Nb 


53 


5.0 


9.0 


7.24 


20 


17 


MoGe 


64 


2.0 


6.25 


160 


2.0 


1.6 


MoGe 


64 


5.0 


6.25 


160 


2.1 


1.6 


MoGe 


64 


7.0 


6.25 


160 


2.0 


1.6 



TABLE I: Transport and superconducting properties of the 
Nb and Moo.yGeo.s samples. Here ds and Ws are the film 
thickness and bridge width respectively, Tc is the sample crit- 
ical temperature, p is the measured specific resistance, JdpiO) 
and Jdp^(0) are extrapolated and calculated critical current 
density at zero temperature. 



III. RESULTS AND DISCUSSION 

To illustrate the raw data, experimentally determined 
values of Jdp as a function of reduced temperature t for 
two bridges of Nb (d^ = 20 nm, w^ = 1 /j,m) and a- 
Moo.yGco.s {dg =64 nm, w^ = 2 /ini) are shown in Fig. 4. 

Between t = 1 and t = 0.85 both curves show a clear 
upturn, which indicates the expected GL behavior. Plot- 

2/3 

ting as a function of t in this temperature region 
results in a straight line, which can be used to extrapo- 
late 3dp{t) to zero temperature. Table 1 shows the val- 
ues of Jdp(O) for all samples investigated. It can also be 
used to obtain the normalized temperature dependence 
(idpit) / JdpiO))'^^"^ , which has a universal form in KL the- 
ory. Plots of this quantity for samples with different 
bridge width are shown in Fig. 5 for Nb and in Fig. 6 
for a-Moo.7Geo.3. Both the absolute values of 3dp{0) and 
the temperature dependence can be directly compared to 
the KL results, which we now briefly reiterate. 
Close to Tc the depairing current density can be written 
as follows : 



Jf/(t) = lMx'/\p)eN{0)vFkBTc{l - T/T^f'^ (1) 

where x{p) is the Gor'kov function controlled by a dimen- 
sionless parameter characterizing the amount of electron 
scattering, p=(hvp)/{2TTkBTc£), with £ the elastic mean 
free path and N(0) the density of states at the Fermi level 
for each spin direction. For £ <C (dirty limit) p oo, 
which yields for xip) ~* 1.33^/^o- Thus, at zero temper- 
ature the extrapolated depairing current density Jdp^{0) 
becomes 



jg'^{0) = 1.26eNiO)vFA{0)^- (2) 

Because of the small mean free path in both types of 
samples, we may assume applicability of the free-electron 
model, so the density of states N(0) can be expressed as 

N{0) = ileWpir' (3) 
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Substituting this formula in Eq. 12 with = hvF/TTA{0) 
and A(0) ~ I.TG/cbTc we obtain 



Jf,^(0) = 244 



1/2 



(4) 



This result is similar to the one obtained in ja] and 
0. Eq. H contains only experimental quantities and 
the pi product, which is known for both materials from 
literature [1,0, Looking now at Figs. 5 and 6, all 



KL- Theory 

Gl Theory 

— •— Nb Vm X 20 nm 
Nb 2nm X 40 nm 
-♦- Nb 2. Sum X 53 nm 
— *— Nb 5|im X 53 nm 




FIG. 5: Experimental results for the pair-braking current den- 
sity J dp normalized to its extrapolated value J dp (0) as a func- 
tion of reduced temperature in Nb bridges of different width 
and thickness as denoted. The black solid and dashed lines 
indicate KL and GL results respectively. 

curves follow GL behavior down to about t = 0.85. 
The values of J(ip(0) extrapolated from this region can 
be compared to the values calculated from Eq. ^ for 
J^^(O). This comparison is made in Table 1 which 
gives all relevant parameters for the different samples. 
Basically, we find quite good agreement for all sample 
widths. In the case of Nb, the most serious deviation is 
found for the 5 ^m bridge, which is presumably due to 
contact heating as a result of the larger current. It is 
interesting to note that the extrapolated values are the 
same as found by Geers et al. who used continuous 
currents and larger bridge widths. The differences are 
in the extent of the GL-regime, which was only found 
down to t = 0.93 in the earlier experiments, and also 
in the temperature dependence below the GL regime. 
There, the temperature dependence is described by the 
full KL-calculation, which was also performed in ref. 0. 
For a single superconducting film, the results for Nb are 
shown in Fig. 5 by the solid line. The smallest sample 
(d=20 nm, w=l pin) follows the KL theoretical curve 
down to t — 0.2 without a significant deviations. Wider 
bridges show a suppression of Jdp{t) with respect to the 
calculated value, again in agreement with earlier results 
0. Presumably, sample heating via contacts and vortex 



fiow occurs even for the short time of a current pulse. It 
appears therefore that using low (pulsed) currents, Jdp(t) 
can be determined correctly over the full temperature 
range for other materials than Al. Circumstances can 
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FIG. 6: Experimental results for the pair-braking current den- 
sity Jdp normalized to extrapolated value Jdp(O) as a func- 
tion of reduced temperature in Moo.yGeo.a bridges of different 
width and thickness as denoted. The black solid and dashed 
lines indicate KL and GL results respectively. 

be somewhat less favorable, however, as shown by the 
measurements on a-Moo.yGeo.s. These were performed 
only for a film thickness of 64 nm. In the GL regime 
the difference between measured and calculated values 
of Jdp{0) is somewhat larger than for Nb (see Table 
1), with the measured values larger than the calculated 
ones. It will be clear that this cannot be due to pile-up 
of current at the samples edges, which would yield the 
opposite effect. Moreover, for amorphous materials this 
should be less of a problem, since the penetration depths 
are very large and actually of the order of the smallest 
bridge width. The difficulty rather lies in the correct 
determination of Jdp(t) close to Tc, with more scatter 
in the individual points. One reason for this may be 
the very low vortex pinning which is characteristic of 
amorphous materials 12, 13|. Another may be that the 
processing of the film during the structuring process 
may lead to changes in the material. For instance, the 
specific resistance we find for the bridges is about 10 % 
lower than for wider structures |l4j . Also, thinner films 
showed increasing p and decreasing Tc, which in this 
thickness regime cannot be well explained by the onset 
of localization effects Q . Since amorphous materials are 
very sensitive to recrystallization, this may be playing a 
role. Still, the difference between Jdp{0) and J^^(O) is 
only 20 %, which may still be considered very good. For 
the temperature dependence (Fig. 6) the result is also 
similar to Nb. For the smallest bridge, the experimental 
curve shows good agreement with the theoretical pre- 
diction, while for wider bridges the values remain too low. 

In summary, we have shown that measurements of 
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depairing currents in conventional type-II superconduc- 
tors with cross-section larger than their characteristic 
lengths and is well possible by using a pulsed 
current method. Using two different superconductors 
with quite different values of their depairing current, 
we found good agreement between experiments and 
theory with respect to both the absolute values and the 
temperature dependence, over essentially the full range 
of temperatures. Such an unambiguous determination 
of a quantity which directly measures the superconduct- 
ing order parameter should also find use in problems 
posed by hybrid systems; in particular, it will be the 
correct quantity to gauge the effects of suppression of 
superconductivity by the injection of (spin-polarized) 



quasiparticles. 
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